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I. INTRODUCTION 

A narrow charmonium-like state X(3872) was ob- 
served in 2003 in the exclusive decay process — >■ 
K^iT+TT-J/tp [H. The A(3872) decays into tt+tt- 7/0 and 
has a mass of mx = 3872.0 ± 0.6(stat) ± 0.5(syst) very 
close to the M^o + Md,o = 3871.81 ± 0.25 mass thresh- 
old [2|. Its width was found to be less than 2.3 MeV 
at 90% confidence level. The state was confirmed in B- 
decays by the BABAR experiment [1] and in pp produc- 
tion by the Tevatron experiments CDF and D0 Q. 
The most precise measurement up to now was done in 
with Mx = 3871.61 ± 0.16 ± 0.19. The new average 
mass given in |3] is 



Mx = 3871.51 ±0.22 MeV. 



(1) 



From the observation of the decay A (3872) J/i'l 
reported by both the Belle ^] and BABAR 8] collabora- 
tions and from the angular analysis performed by CDF 
0, it was shown that the only quantum numbers com- 
patible with the data are J^'~^ — 1++ or 2 '". How- 
ever, the observation of the decays into D^D 7r° by the 
Belle and BABAR collaborations allows one to 

exclude the choice 2 ^ because the near-threshold decay 

tt" is expected to be strongly suppressed for 

J = 2. 

The Belle collaboration has reported evidence for the 
decay mode X — > tt+tt^tt^ J/t/; with a strong three-pion 
peak between 750 MeV and the kinematic limit of 775 
MeV |7|] , suggesting that the process is dominated by the 
sub-threshold decay X — >• wJ/^. It was found that the 
branching ratio of this mode is almost the same as that 
of the mode X — )■ 7r~^Tr~ J/ip: 



B(X J/7/;7r+7r-7r") 
B{X J/7/;7r+7r-) 



1.0±0.4(stat)±0.3(syst). (2) 



ate w-meson with isospin whereas the two-pion decay 
proceeds via the intermediate p- meson with isospin 1. 
Also the two-pion decay via the intermediate p-meson is 
very difficult to explain by using an interpretation of the 
A (3872) as a simple cc charmonium state with isospin 0. 
There are several different interpretations of the A(3872) 
in the literature: 

• The A(3872) is a molecule bound state (D°D*°) 
with small binding energy (see, Refs. ([l2|-[25j)): 

• The A (3872) is a tetraquark state composed of a 
diquark and antidiquark (see, Refs. ([26i]-[31i]')'): 

• The A (3872) should be interpreted in terms of 
threshold cusps [32[; 



These observations imply strong isospin violation be- 
cause the three-pion decay proceeds via an intermedi- 



• The A (3872) should be interpreted in terms of hy- 
brids ^33] and glueballs (ssj . 

A description of the current theoretical and experimental 
situation for the new charmonium states may be found 
in the reviews [36j-[43|. 

The most intriguing question at present is whether the 
X(3872) is a losely-bound charm-meson molecule with a 
binding energy of Mx - (M^.o + Md<^) = -0.30 ± 0.40 
MeV, or a tetraquark composed of a color diquark and 
color antidiquark. Recently, the molecular nature of the 
the X(3872) has been carefully discussed in the literature. 
The authors of [H] assumed that the X(3872) is a D^D* ° 
molecule and, based on this assumption, estimated its 
prompt production cross section at the Fermilab Teva- 
tron. They presented a theoretical upper limit on the 
prompt production cross section which is about 30 times 
smaller than the observed prompt production rate at 
the Tevatron. The conclusion was that S-wave resonant 
scattering is unlikely to allow the formation of a loosely 
bound D'^D* " molecule in high energy hadron collisions. 
In reference |46j| it was argued that one needs to take into 
account charm- meson rescattering in the analysis of [45| . 
As a consequence the theoretical upper limit is increased 



2 



by orders of magnitude and is then compatible with the 
observed production rate at the Tevatron. This conclu- 
sion was later on criticized in [13] ■ First, using the results 
of [46*1, the authors of argue that a new unobserved 
Xs(4080) molecule composed of a DsD* pair should have 
been observed at the Tevatron. Second, they cast some 
doubts on the applicability of the Watson theorem for 
final state interactions in the above calculation. 

The tetraquark state interpretation of the X(3872) 
was successfully applied to describe the available exper- 
imental data on the decays X — >■ J/^tt+tt" and X — > 
J/07r+7r~7r'' [2^. By using an effective three- meson La- 
grangian with the coupling taken from a similar analysis 
of light scalar mesons, the authors found the width of the 
X- meson to be 1.6 MeV in accordance with the experi- 
mental bound. Contrary to this, the values of the widths 
calculated by using QCD sum rules were found to be too 
large - around 50 MeV [2^ . 

In this paper we provide an independent analysis of 
the the properties of the X(3872) meson which we inter- 
pret as a tetraquark state as in [2^. We work in the 
framework of the relativistic constituent quark model 
which has recently been extended to include infrared 
confinement effects [1^. The improved model (481 is 
a successful generalization of the relativistic constituent 
quark model which some of us have developed over many 
years [l^-js^. The relativistic constituent quark model 
can be viewed as an effective quantum field theory ap- 
proach to hadronic interactions based on an interaction 
Lagrangian of hadrons interacting with their constituent 
quarks. Once the relevant interpolating quark current 
is written down one can evaluate the matrix elements 
of the physical processes in a self-consistent way. The 
nice feature of this approach is that multiquark sys- 
tems as e.g. baryons and tetraquarks can be treated on 
the same footing as the simplest quark-antiquark states. 
The coupling strength of hadrons with their interpolating 
quark currents is determined by the compositeness con- 
dition Zh = [sij where Zh is the wave function renor- 



malization constant of the hadron. Matrix elements are 
generated by a set of quark loop diagrams according to a 
1 /Nc expansion. The ultraviolet divergences of the quark 
loops are regularized by including vertex form factors for 
the hadron-quark vertices which, in addition, describe 
finite size effects due to the non-pointlike structure of 
hadrons. The relativistic constituent quark model con- 
tains only a few model parameters: the light and heavy 
constituent quark masses, the confinement scale, and the 
size parameters that describe the size of the distribution 
of the constituent quarks inside the hadron. 

The paper is organized as follows. In Sec. II we con- 
struct the nonlocal generalization of the four-quark in- 
terpolating current of the X(3872) written down in [26j |. 
This leads to a nonlocal effective Lagrangian describing 
the interaction of the X(3872) meson with its constituent 
quarks. The coupling strength of the X(3872) w.r.t. its 
constituent quarks is determined from the compositeness 
condition Zh = 0. We also briefly discuss the imple- 
mentation of infrared confinement in our scheme. In 
Sec. Ill we calculate the matrix elements of the tran- 
sitions X ^ J/ip + p{w) and X D + D* . The re- 
sults are then used to evaluate the widths of the decays 
X J/iP + 2Tr{3Tr) and X ^ D° + D° + Tr°. In Sec. IV we 
present the results of our numerical analysis and compare 
our results with the results of other approaches. Finally, 
in Sec. V we summarize our findings. 



II. THEORETICAL FRAMEWORK 
A. Effective Lagrangians 

The authors of [1^ suggested to consider the X(3872) 
meson as a J^'-^ = 1+"'" tetraquark state with a sym- 
metric spin distribution: [cq]s=o [eg] 5=1 -f [eg] 5=1 [cq]s=o, 
{q = u,d). The nonlocal version of the four-quark inter- 
polating current reads 



Wi 



W2 = Wc 



dxid ix -^^WiXij ^x(^'^{xi - Xj)"^^ (3) 

V i=l ) i<j 

I [qaiXi)C-f^Cbixi)][qd{x3h^CCeix2)] + [qa{x4)C-f'' Cbixi)][qdix3)-f^ CCeix2)] }, 

Too 



nic 



2{mq + rUc) ' 



2{mq -I- TOc) ' 



r 



Feynman diagrams. The effective interaction Lagrangian 
describing the coupling of the meson Xq to its constituent 



The matrix C = 7*^7^ is the charge conjugation ma- 
trix: C = = c-i = -C^, CT^C-i = ±r, (" + " 
for V = S,P,A and " — " for F — V, T). The numbering quarks is written in the form 
of the coordinates Xi is chosen such that one has a con- 
venient arrangement of vertices and propagators in the £. 



gxXq^{x) ■ (x), {q = u,d). (4) 
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The state X„ breaks isospin symmetry maximally: 



V2\ ^/2 + V2 



(5) 



7=0 



1=1 



The authors of '26'] take the physical states to be a linear 
superposition of the Xu and Xd states according to 

Xi = Xiow = Xu cos 9 + Xd sin 9, 
Xh = Xhigh = -Xu sai9 + Xd cos 9. (6) 

The mixing angle 9 can be determined from fitting the 
ratio of branching ratios Eq. ([2|) . 

The coupling constant gx in Eq. (|4]) will be deter- 
mined from the com pos iteness condition Zh = 0, see, 
e.g. Refs. [131, IHll' nil- The compositeness condition 
requires that the renormalization constant Zx of the el- 
ementary meson X is set to zero, i.e. 



Zx = l-U'ximj,)=0, 



(7) 



where Ilx{p^) is the scalar part of the vector- meson mass 
operator 



1 / Pfj^Pi. 

9^cly — 



(8) 



The Fourier transform of the vertex function 
j)^) can be calculated by using appropi- 



ately chosen Jacobi coordinates. One has 

W3 - W4 



2W2 + W3 + W4 

2wi + W3 + W4 

X2 = X ^1—^ Pi 



2V2 



Wi — W2 
X3 = X Pi 



2V2 



Wi - W2 

X4 = X — Pi 



2V2 

W3 - W4 
2V2 

Wi + W2 + 2W4 

2V2 

Wi + W2 + 2W3 



2V2 



2V2 



P2 + 
P2 + 
P2 - 
P2 - 





+ 






2 




W3 


+ 


W4 




2 




Wi 


+ 


W2 




2 




Wi 


+ 


W2 


2 



P3, 
P3, 
P3, 
P3, 



where x — XiWi and ^ (a 



then has 



i<i<j<i 



,)2=EP?-0ne 



x\ 7 Xx%-Xj) ^ 

X ^xiPl 



{2^Y 



,P4) 



\i=i 



mjdp. 



i Y2 UJiPi „ 

3 ^x{R'), 



(9) 



3 3 
where il^ = J2 '^f ^-^i^ ^ J2 Pi- The Jacobi coordi- 

nates in momentum space read 



Pi ~P2 
2V2 ' 

Pi +P2 



UJ2 



P1+P2 + 2p3 

2V2 



(10) 



For calculational convenience we will choose a simple 
Gaussian form for the vertex function <^x{—^^)- The 
minus sign in the argument of this function is chosen to 
emphasize that we are working in Minkowski space. One 
has 



exp {nyAj,) 



(11) 



where the parameter Ax characterizes the size of the X- 
meson. Since fi^ turns into — Jl^ in Euclidean space the 
form pTjl has the appropriate fall-off behavior in the Eu- 
clidean region. We emphasize that any choice for $x 
is appropriate as long as it falls off sufficiently fast in 
the ultraviolet region of Euclidean space to render the 
corresponding Feynman diagrams ultraviolet finite. As 
mentioned before we shall choose a Gaussian form for 
<f>x for calculational convenience. 

We are now in the position to write down an explicit 
expression for the derivative of the mass operator appear- 
ing in Eq. ([7]). The corresponding three- loop diagram de- 
scribing the X-meson mass operator is shown in Fig. [T] 





FIG. 1: Diagram describing the X„-meson mass operator. 
One has 



2 3i 



3p2 \9^- „2 



n 



\J {2^^ 



2l„,5 



tr 



X I — Wctr 
— Wctr 

+«;,tr [5p]^55[2]^5] [5[3]^M5[13] } ( 

where we have introduced the abbreviations 



12) 
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S^^l = S',(A;2 + Wqp), 



In the next section we shall describe how to evaluate the 
integral P^ . 



a quadratic form of the external momenta. Tensor loop 
integrals are calculated with the help of the differential 
representation 



1 d 



2di 



2kr 



(16) 



2 ^ 



We have written a FORM [61| program that achieves the 
necessary commutations of the differential operators in a 
very efficient way. After doing the loop integrations one 
obtains 



00 



,/3n): 



(17) 



B. Infrared confinement 

In [4^ we described how to integrate n-point one-loop 
diagrams and how to implement infrared confinement of 
quarks in this process. Since the present application in- 
volves also multi-loop diagrams we need to extend our 
loop integration techniques to the case of arbitrary num- 
ber of loops. Let n, £ and m be the number of the prop- 
agators, loops and vertices, respectively. In Minkowski 
space the ^-loop diagram will be represented as 

n{pi, ...,pm) = 



11=1 



^2 _ -^7/'^^ ^2 



(13) 



where the vectors ki are linear combinations of the loop 
momenta ki. The Vi are linear combinations of the ex- 
ternal momenta pi to be specified in the following. The 
strings of Dirac matrices appearing in the calculation 
need not concern us since they do not depend on the 
momenta. The external momenta pi are all chosen to be 

m 

ingoing such that one has Pi =0. 

Using the Schwinger representation of the local quark 
propagator one has 



S{k) = (m+ /^) yd/3e-^ 




(14) 



For the vertex functions one takes the Gaussian form. 
One has 



=exp i = l,...,m, (15) 



where the parameters /3i+„ — s,; — 1/A| are related to 
the size parameters. The integrand in Eq. (jl3p has a 
Gaussian form with the exponential kak + 2kr + R where 
a is £ X ^ matrix depending on the parameter /J^, r is the 
f-vector composed from the external momenta, and R is 



where F stands for the whole structure of a given dia- 
gram. The set of Schwinger parameters Pi can be turned 
into a simplex by introducing an additional t-integration 
via the identity 



1 = J dt5{t-Y,P^) 



(18) 



leading to 



n= dtt 



d"a(5(l-^a, 



F(toi,...,to„). (19) 



There are altogether n numerical integrations: {n — 1) 
Qf-parameter integrations and the integration over the 
scale parameter t. The very large i-region corresponds 
to the region where the singularities of the diagram with 
its local quark propagators start appearing. However, 
as described in [48| , if one introduces an infrared cut-off 
on the upper limit of the t-integration, all singularities 
vanish because the integral is now convergent for any 
value of the set of kinematic variables. We cut off the 
upper integration at 1/A^ and obtain 



l/A^ 1 

= j dte~^J d"aS(l - J2 Pi*' 



ai,...,to„). (20) 



By introducing the infrared cut-off one has removed all 
potential thresholds in the quark loop diagram, i.e. the 
quarks are never on-shell and are thus effectively con- 
fined. We take the cut-off parameter A to be the same 
in all physical processes. The numerical evaluations have 
been done by a numerical program written in the FOR- 
TRAN code. 

As a further illustration of the infrared confinement 
effect relevant to the applications in this paper we con- 
sider the case of a scalar one-loop two-point function. 
One has 



^p' 



-ski 



5 



where we have collected all the nonlocal Gaussian vertex 
form factors in the numerator factor 6"**^^=. Note that 
the momenta ks, Pe are Euclidean momenta. Doing the 
loop integration one obtains 



U^ip') = dt 



t 



zioc = TO^ - a(l - 



da exp 



- tz 



Zl 



loc 



Zl 



st 
s + t 

a-i)V.(21) 



The integral 112 (p^) can be seen to have a branch point 
at p^ = 4m^ because zioc is zero when a — 1/2. By 
introducing a cut-off on the ^-integration one obtains 



dt 





da exp 



loc 



st 

s + t 



Zl 



(22) 



The one- loop two-point function Il2{p^) Eq. (|22|) can be 
seen to have no branch point at p^ = Arn? . 



III. THE TRANSITIONS X_^ J/ip + p(lj) AND 
X ^ D + D' 

In this section we evaluate the matrix elements of the 
transitions X ^ J/tp + p{uj) and X D + D* . The 
relevant Feynman diagrams are shown in Fig. [5] Since 




FIG. 2: Feynman diagrams describing the decays X — >■ J/'^-\- 
p{uj) and X ^ D + D*. 

the X(3872) is very close to the respective thresholds in 
both cases, cif. 



mx - (mj/tl) + nip) 



-0.90 ± 0.41 MeV, 
-0.30 ± 0.34 MeV 



the intermediate p, uj and D* mesons have to be treated 
as off-shell particles. 



(x„ (p, A* ) J/V^ (gi , i^) -f t;" ((72 , P)) 



6 9x 9j/^ 9v° 



d*ki 



(27r)4i J (27r)4i 
(fci-}-iqi)2)$„o(-(fc2 + ig2) 



X tr 



ij^SMYSciki + qi)rSu{k2)YSu{k2 + <Z2) 



(23) 



Kf = |(fci + igi)2 + i(/c2 + i(?2)^ + 3(w„gi-Wcg2)^- 

where — p,uj. In the case where the X and 7/0 are on 
mass-shell, i.e. C/i (91+92) — ^^^^ e^^gi = 0, the number 
of independent Lorentz structures reduces to 6. Due to 
the obvious relations: 



-M{Xu^ m + p), 



M{Xd ^ J/i^ + p) 
M{Xd Jli) + uj) 

one can express the decay amplitudes of the physical 
states defined by Eq. ([6]) via the decay amplitudes of 
the Xu as 

M{Xi Jli) -I- w) = (cos B + sin d) M{X^ J/i/; -I- w) , 
M{Xh J/i^ + w) = (cos 6 - sin 9) Af (X„ ^ J/i/> + u) , 
M(X, -> J/V^ + p) = (cos e - sin 9) M{Xu ^ J/i/- + p) , 
M(X,, ^ J/t/- + p) = -{cos 9 + sine) M(X„ ^ J/^j + p). 

Next we turn to the decay X — > l) + D* . The generic 
matrix element is written as 

D{qi) + D*{q2,y) 



M^^(x,{p,p) 
= 3V29X 9d 9d* 



d^k^ f d^k2 



X tr 



j {2n)H^''\ 

5 ( - (fci -I- Wc^l)^) ( - (^2 + Wc92)^) 

755,(fci)7^5,(fci -I- 9i)7'^^c(fe)7''^,(fc2 + 92) 



-I- (rriq O TTic, Wq O Wc) 



r(i) 



(2) 



r(3) 



9 ^^^XDD- + 9l 91 ^^^XDD' + 9l 92 ^^^XDD 



+ 92 91 '^^XDD' + 92 92 ^^^XDD* ■ 
Kl = \{ki - k2f + \{ki - k2 + qi - q2f 

+ \{ki+k2+Wcpf. 



(24) 
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Note that the decay of the X-meson into DD is forbid- 
den. 



IV. NUMERICAL ANALYSIS 



A. X-decays 



Using the matrix elements ([25]) for the decay X — > 
J/ip + p{uj) one can evaluate the decay widths X — > 
J/t/j + 27r(37r). We employ the narrow width approxi- 
mation which was extensively discussed in Ref. [5a] and 
also used in Ref. [59|. One has 



Here p*{(f') = \^/'^{'m\,rn?ji^,q'^) /2mx is the momen- 
tum of the J/Tp or the (n7r) system in the center-of-mass 
frame and (nm,r)^ < <f < {mx—fnj/^Y defines the kine- 
matic region of the respective processes with n = 2 for the 
p meson and n = 3 for the lo meson. The essential point 
which allows one to derive Eq. ((25)) is the omission of 
polarization correlations. As was shown in Ref. [55l | they 
produce no effects if the intermediate state is on-shell. 
Note that in our calculation we keep the q^— dependence 
of the matrix elements as given by Eq. (l23l) . We use 
the masses and widths of the p(a;)-mesons from 0] (all 



in MeV): 



782.65, 



= 775.49, 
= 8.49. 



146.2, B{p - 
Stt) 0.892, 



27r) 



dT{X ^ J/iIj + nTT) 



dq2 



^ ■\\Mxjv\'' 



8m^7r 3 



\M 



XJv 



pol 



^ nTT), (25) 



The adjustable parameters of our model are the con- 
stituent quark masses mq , the scale parameter A charac- 
terizing the infrared confinement and the size parameters 
Am- They were determined by using a least square fit to 
a number of physical observables, see 48]. Below we dis- 
play the numerical values for the parameters which are 
relevant to the present paper. 





rric A 


A^ 


Ap/^ 


Ad 


Ad- 


Aj/^ A^^ 




0.217 0.360 


1.6 0.181 


0.711 


0.295 


1.4 


2.3 


3.3 3.0 


GeV 



(26) 



Note that our fit values for the size parameters of the 
pion and the p-meson are in qualitative agreement with 
those found in the quark model based on the instanton 
vacuum Ref. (57| . 

There are two new free parameters: the mixing angle 
9 in Eq. (|6]) and the size parameter Ax- We have varied 
the parameter Ax in a large interval and found that the 
ratio 



T{Xu J/i/' + Stt) 
r(Xu J/tp + 2Tr) 



0.25 



(27) 



is very stable under variations of Ax- Hence, by using 
this result and the central value of the experimental data 
given in Eq. ([5]), one finds 



TjXiM J/^ + 3 7t) 
r{Xi,h^ J/i: + 2TT) 



0.25 



(i 



1 ± tan 



q= tan0 



1 (28) 



which gives « ±18.4° for Xi {"+") and X,, ("-"), re- 
spectively. This is in agreement with the results obtained 
in both 6* w ±20° and [29]: 6 w ±23.5°. The decay 
width is quite sensitive to the change of the size parame- 
ter Ax - A natural choice is to take a value close to Aj/^ 
and A^^ which are both around 3 GeV. We have varied 
the size parameter A^ from 3 up to 4 GeV and display the 
dependence of the decay width in Fig. [3] One can see that 
the decay width T{X — >■ J/0 + n tt) decreases from 0.30 up 
to 0.07 MeV, monotonously. This result is in accordance 



with the experimental bound r(X(3872)) < 2.3 MeV and 
the result obtained in [2^: 1.6 MeV. 

In a similar way we calculate the width of the decay 
X — > D^D^TT^ which was observed by the Belle Coll. and 
reported in 10]. Again using the narrow width approxi- 
mation the differential rate reads 



dr(X„ ^ i: >"£)"7r"^ 



1 1, 



2ni\TT 3 



IM. 



XDD* 



^ r^,om^,o p*(g2)g(i^*"-» J"^°) ^29) 



IM, 



XDD* 



_„2\2_|_p2 2 



\ " ipM ^i- i\,r |2 



3 i^x 

pol 



where the matrix element is defined by Eq. ([24]) 

andp*(<?2) = A^/2(m^,TO^o,9^)/2mx is the momentum 
in the center-of-mass system. The kinematical region de- 
fined by 

(mjjo + m^o)^ < 9^ < i^x — fTT-Do)^ 

is very narrow 3.99928 < < 4.02672 GeV^ Note that 
we have taken into account both channels with the inter- 
mediate and D*'^ mesons. We use the masses and 
widths of the D* + and D* " mesons given in 0] and [l6| . 
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[27| (all dimensional quantities in MeV): 



2010.27, 
2006.97, 
0.070, 



Td-+ = 0.096, 
B{D*+ ^ D+TT° 
B{D*" 



0.307, 
0.619. 



Keeping in mind that 

T{Xi D°D\°) = cos6l2r(X„ ^ D°D\°), 

we have varied Ax from 3 up to 4 GeV and found that 
the decay width T{Xi -> Z)°£)°7r°) decreases from 1.88 
up to 0.41 MeV, monotonously. In Fig. [3] we plot the 
dependence of the decay widths T{Xi — >■ D^D^tt'^) and 
r(X — >■ J/'it+ rnr) on the size parameter Ax- Using the 
results of [2] 

10^S(B± ^ if±A:) • B{X J/^tt+t:-) 0.95 ± 0.19, 
10^6(5=^ -> isT^X) • B{X D°D"tt°) = 10.0 ± 4.0 



one calculates the rate ratio 

r{X D°D^TT°) 



r(A: -> j/i/^TT+TT-) 



= 10.5 ±4.7 



(30) 



The theoretical value for this rate ratio depends on the 
size parameter Ax as Fig. |3] shows. One has 



T{X D^D^TT^-^ 



T{X J/iPtt+tt-) 



6.0 ±0.2. 



thcor 



(31) 



where the theoretical error reflects the Ax dependence of 
the ratio. The ratio lies within the experimental uncer- 
taintities given by Eq. (pO| . 




FIG. 4: Diagram describing the X-resonance contribution to 
the J/(/)-dissociation process. 

Charged D-mesons 




Neutral D-mesons 




FIG. 3: The dependence of the decay widths T{Xi 
D'^ D°-K^'') and r(X J/%1) + nvr) on the size parameter Aji 



B. J/(/;-dissociation 

The last topic which we would like to discuss is the im- 
pact of the intermediate X-resonance on the value of the 



FIG. 5: The cross sections of the processes J/tp -I- — !> X — > 
D + D* . Charged D-mesons- upper panel, neutral D-mesons- 
lower panel. 



J/0-dissociation cross section, see [5^-[60|. The relevant 
s-channel diagram is shown in Fig. |4l 

We evaluate the J/V'-dissociation cross sections without 
using the narrow-width approximation. One has 



D + D* 



= 2{coseTsme)^a{J/i} + v° X^ ^ D + D*), 



a{J/ip + v° Xu-^ D + D*) 

16 TT s AV2(s, m2/^, ™^o) 9 ^ (s - mj,)^ + r|m| ' 
A = e^ji^el, AV,( - 5^" + t^)e%.M^, (32) 
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where p = pi + P2 = 9i + 92- = p (minus sign) 
or uj (plus sign). We also neglect isotopic spin break- 
ing effects. Note that E = > 111^+ + TO£).+ for 
charged D- mesons in the final states and E > mj/^ + m-^o 
for neutral D-mesons in the final states. In the first 
case the cross section is zero at the threshold of reac- 
tion E = m£)+ + mjjt + . In the last case the cross sec- 
tion blows up at £^ = mj/^ + rriyo because the channel 
J/tp -I- — > + D* ° is exothermic and the kinematical 
function A^/^(s, mjy^^, m^o) in the denominator is equal 
to zero at this point. We take Tx = 1 MeV in the Breit- 
Wigner propagator and set Ax =3.5 GeV when calcu- 
lating the matrix elements. We plot the behavior of the 
relevant cross sections in Fig. [S] One can see that in 
the case of charged D-mesons (upper panel in Fig. [5]) the 
maximum value of the cross section is about 0.32 mb at 
E = 3.88 GeV. This result should be compared with the 
result of the cross section a{J/ip + TT — )• D + D*) w 0.9 mb 
ai E = 4.0 GeV, see, [5S| and the result of the cross sec- 
tion a{J/ii} + p D + D*) ^2.9 mh at E = 3.9 GeV, see, 
[ssj . Thus the X- resonance gives a sizable contribution 



to the J/0-dissociation cross section. It would be inter- 
esting to do a complete analysis of the J/ip dissociation 
cross section in view of our new results on the s-channel 
contribution of the X(3872) tetraquark state. However, 
this certainly is beyond the scope of the present paper. 
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